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ABSTRACT 

We investigate from the fractal viewpoint the way in which the dark matter is grouped 
at z — in the Millennium dark matter cosmological simulation. The determination of the 
cross to homogeneity in the Millennium Simulation data is described from the behaviour of 
the fractal dimension and the lacunarity. We use the sliding window technique to calculate 
the fractal mass-radius dimension, the pre-factor F and the lacunarity of this fractal relation. 
Besides, we determinate the multi-fractal dimension and the lacunarity spectrum, including 
their dependence with radial distance. This calculations show a radial distance dependency of 
all the fractal quantities, with heterogeneity clustering of dark matter haloes up to depths of 
100 Mpc/h. The dark matter haloes clustering in the Millennium Simulation shows a radial 
distance dependency, with two regions clearly defined. The lacunarity spectrum for values of 
the structure parameter q ^ 1 shows regions with relative maxima, revealing the formation 
of clusters and voids in the dark matter haloes distribution. With the use of the multi-fractal 
dimension and the lacunarity spectrum, the transition to homogeneity at depths between 100 
Mpc/h and 120 Mpc/h for the Millennium Simulation dark matter haloes is detected. 

Key words: (Cosmology:) dark matter - large-scale structure of the Universe - methods: 
statistical 



1 INTRODUCTION 

The current paradigm, which is the basis of the standard cosmo- 
logical model, states that the large-scale universe is homogeneous 
and isotropic, i.e., there is no privileged place or direction in the 
universe (Martinez & Saar 2002) in a way that the observed inho- 
mogeneities are local in character and they should vanish at suffi- 
ciently large scales (Peacock 1999; Longair 2008). This principle 
is based mainly on the observations of the cosmic microwave radi- 
ation whose isotropy is one part in a hundred thousand (Baryshev 
& Teerikorpi 2005), and philosophical considerations: the observed 
universe should be the same for any observer, the Copernican Prin- 
ciple. 

Despite the success of the physical models constructed from 
this cosmological principle, some researchers now are debating 
heavily whether or not the observations from more recent galaxy 
catalogs corroborate the assumption of homogeneity and isotropy 
at scales greater than 60 Mpc/h (Sarkar et al. 2009; Capozziello 
& Funkhouser 2009), or if instead the galaxies are grouped into 
highly structured hierarchical pattern, with properties of fractality 
(Gabrielli et al. 2005; Verevkin, Bukhmastova & Baryshev 201 1). 

In the last decade, the interest in the fractal analysis of the 
galaxy distribution has yielded promising results, becoming a fer- 
tile field for cosmological physics research (Baryshev & Teerikorpi 



2005). In particular, we do not know if the cold dark matter dis- 
tribution have a fractal behaviour in the same way that is observed 
in the galaxy clustering. Some authors propose the same fractal di- 
mension for the galaxy distribution and for the dark matter, where 
galaxies act as tracers of dark matter (Baryshev 2008), while other 
researchers assign different fractal dimension; D=3 for dark matter, 
i.e. homogeneity in accord with cosmological principle and D=2 for 
galaxy clustering (Durrer & Sylos Labini 1998). Authors like Gaite 
(2005) suggest the use of multi-fractal models in order to analyse 
the dark matter haloes in numerical N-body simulations. He found 
that haloes of similar mass have a fractal distribution with a given 
dimension that grows as the mass diminishes. 

The homogeneous self-gravity matter can develop fractal 
structure (Combes 1998). For this reason the N-body cosmolog- 
ical simulations can be used like a numerical laboratory where the 
fractal analysis of matter under gravitational interaction can be per- 
formed. Besides, in simulations with sufficient number of particles 
it is possible to detect the fractal behaviour with the scale, includ- 
ing the transition to homogeneity. The statistical algorithms used in 
the analyses of the simulation data can be applied subsequently to 
galaxies surveys (Gabrielli et al. 2005). 

In this paper we present the calculations of mass-radius fractal 
dimension, mono-fractal lacunarity, multi-fractality and lacunarity 
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spectrum of dark matter haloes distribution from the Millennium 
Simulation Project. We propose the use of multi-fractal dimension 
with the complementary measure of the lacunarity spectrum as a 
method for determining the homogeneity scale, i.e., the distance 
beyond which the transition to homogeneous regime should occur 
in a fractal matter clustering. 

The paper is divided as follows: Section 2 presents the mathe- 
matical foundations involved in the fractal approach to the cluster- 
ing of large scale matter; section 3 summarises the characteristics of 
the Millennium simulation, focusing on the cosmological parame- 
ters and scale used as the basis of the simulation; section 4 discusses 
the calculation of mass-radius fractal dimension and lacunarity as 
functions of radial distance in the clustering of dark matter haloes 
from the Millennium Simulation; section 5 deals with the calcula- 
tion of the multi-fractal spectrum and the spectrum of lacunarity 
on the same data from the Millennium Simulation; in Section 6 the 
discussion of the fractal results is presented and finally the conclu- 
sions are given in section 7. 



2 THEORETICAL DEVELOPMENT 
2.1 Dimension Concept 

The most general definition for the dimension measure was devel- 
oped by Felix Hausdorff in 1918 (Hausdorff 1918), looking to de- 
fine the dimension concept for any metric space. He first consid- 
ered the number N(e) of balls of radius e needed to cover a subset 
A completely. The d-dimensional measure is proportional to the 
limit: 



h d (A) 



lim N(e)e a 

e->0 



(1) 



Hausdorff proved that there is only one value d for which this 
is different from zero or infinity. For this value is satisfied that: 



h d (A) = 







if d<D H {A) 
if d> D„{A) 



(2) 



where Dh is by definition the Hausdorff dimension for the subset 
A, with Dh the value for which N(e) oc 1 /e d . In order to decide 
whether a set is fractal or not it is necessary to compare the above 
definition with the topological dimension, defined in a simple way 
as the number of independent directions in which one can move 
around a given point of the set. In this form Mandelbrot define a 
fractal as a set for which the Hausdorff dimension strictly exceeds 
its topological dimension. (Mandelbrot 1983) 

It is not easy to determine the Hausdorff dimension for galaxy 
clustering, although there are methods developed to apply the 
Hausdorff dimension, like the Minimal Spanning Tree (Martinez 
& Saar 2002). For this reason it is necessary to use alternative con- 
cepts of fractal dimension for the analysis of the matter distribution 
on a large scale in the universe. One of the most used concepts is 
the mass-radius fractal dimension. According to Blumenfeld and 
Mandelbrot (Blumenfeld & Mandelbrot 1997), given a sphere of 
radius r in a space with Euclidean dimension d, which encloses a 
self-similar fractal structure, the total mass measure M(r) enclosed 
by the sphere take the form: 



M(r) = Fr L 



(3) 



here M(r) represents the number of galaxies (or in our case dark 
matter haloes) in a sphere of radius r, D m is the fractal mass-radius 



dimension, and F is a constant related to the average distance be- 
tween nearest neighbours. The exponent D m , the fractal dimension 
is smaller than the embedding Euclidean dimension d for inhomo- 
geneity condition in the set. 

Although the mass-radius dimension has the advantage of ex- 
tending the measure of fractal dimension of the set on large scales, 
its accuracy is limited as it does not move the centre of observation, 
skewing the statistical analysis. Therefore it is necessary to refine 
the determination of fractality including average looking at several 
points in the sample as centres. According to Gabrielli et al. (2005) 
it is possible to average the mass-radius fractal dimension moving 
on different occupied points of the galaxy distribution. The average 
mass-radius relation is a fractal measure which reflects the condi- 
tional cosmological principle: The cosmos appears statistically the 
same to all observers situated on a galaxy (point of a fractal) but 
not in a region of void (Mandelbrot 1983). 

To fully describe a fractal it is necessary to specify how the 
pre-factor F behaves. For this reason it is necessary to use a com- 
plementary fractal quantity, the lacunarity, which is a general mea- 
sure for the analysis of spatial patterns. 



2.2 Lacunarity Definition 

The majority of the studies aiming to determine the fractality of 
large-scale structures in the universe begin with the calculation of 
the fractal dimension, but fractals with similar dimensions can oc- 
cupy the space where they are embedded in different ways. There- 
fore, it is necessary to recognise that the fractal dimension is not 
enough to uniquely characterise fractal sets, and for this reason Blu- 
menfeld & Mandelbrot (1997) proposed a complementary fractal 
measure, the lacunarity. 

The lacunarity appears in the context of the mass-radius fractal 
dimension, although we can also extend this concept for other frac- 
tal dimension definitions like the correlation dimension (Martinez 
& Saar 2002). Lacunarity high values suggest the presence of large 
empty regions inside the clustering and therefore an increased het- 
erogeneity, while fractals with low values indicate approximation 
to homogeneity. 

The starting point for calculating the lacunarity is the pre- 
factor F from the mass-radius relation, Eq. (3). There is not a 
unique definition to describe the properties of the factor F, but 
Mandelbrot proposed a description based in series of factors of 
variability (Blumenfeld & Mandelbrot 1997), 



C k (D) 

mD)r 



(4) 



where C k (D) is the kth cumulant, a statistical measure of the cor- 
relation designed to go to zero whenever any one or more quanti- 
ties under study become statistically independent. The cumulants 
are related with the moments, in particular k\ = /ii is the mean, 
ki is the variance, and fc 3 = ((A — /ii)' 1 ) for a variable X under 
study. The simplest of this variability factors, is the second order 
variability factor S2 through which the lacunarity <E> is defined as: 



((F-(F)f) _<F 2 ) 



1. 



(5) 



The Lacunarity analysis complements the studies of fractal 
sets enabling us to distinguish between similar fractal patterns 
which occupy the space in a different form. Martinez & Saar (2002) 
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applied the above definition of lacunarity to the Las Campanas sur- 
vey data and the corresponding fractal Levy flight model. They 
found great difference in lacunarity behaviour with similar fractal 
mass-radius dimension. 

The average mass-radius dimension and the lacunarity are suf- 
ficient to characterise sets where the fractal calculation of matter 
clustering at different scales and levels of mass density are the 
same. In the case of scale dependency, we must go beyond in the 
characterising of the large scale matter distribution. This kind of 
analysis can be performed by the multi-fractal formalism. 



2.3 Generalised Dimension 

The dimension definitions presented above, represent particular 
cases of the multi-fractal spectrum of generalised dimension, 
method applied by the first time to the large scale structure in the 
universe by Coleman & Pietronero (1992). The large scale matter 
distribution can be characterised with a fractal spectrum defined 
from a generalised correlation. First, is necessary to introduce the 
correlation integral C2, a function capable to measure the number 
of neighbours that on average a chosen centre has within a distance 
r. We outline the development in accordance with the notation de- 
scribed in Bagla, Yadav & Seshadri (2008): 



corresponds to an arbitrary real number. From this generalised cor- 
relation it is possible to do an expansion in powers of log(r) as 
described by Provenzale, Spiegel & Thieberger (1997) and to cal- 
culate directly the multi-fractal dimension and the lacunarity spec- 
trum: 



log [C^r) 1 '*'- 1 '] = D q log(r) + log(F q ) + O 



1 

logir) 



(11) 

Keeping only the first two terms on the right side, we have the re- 
lation between the generalised correlation integral and the gener- 
alised fractal dimension: 



C q (r 



,1/(9-1) 



(12) 



In this manner the generalised dimension and the generalised 
lacunarity can be defined in the same way as for the mass-radius 
fractal dimension, the generalised dimension D q : 



D 



1 



dlogC 9 (r) 



(13) 



(q — 1) dlogr 

and the corresponding generalised lacunarity from the pre-factor 
F q for every structure parameter q: 



C„(r) 



1 

~NM 



(6) 



where N is the total number of particles inside the distribution, M 
is the number of particles used as centres and the summation is 
performed over the set of chosen centres. n»(r) is the number of 
particles within a radial distance r from a particle at the point i, 
defined as: 



(7) 



where the summation is performed over all the particles in the sam- 
ple. The coordinates of each particle in our space of three dimen- 
sions are denoted by Xj, and is the Heaviside function, defined 
such that O(x) = for x < and Q(x) = 1 for x ^ 0. The num- 
ber of particles around each centre rii(r) is determined by counting 
the number of particles around that centre that lie inside a comov- 
ing sphere of radius r from it. 

From Equation (6) the correlation dimension is defined simi- 
larly to the definition of mass-radius dimension: 



C 2 (r) 



„°2 



(8) 



So that the correlation dimension is calculated as a derivative: 



D 2 = 



d log C 2 (r) 



(9) 



dlogr 

From the correlation integral C2<Y), Equation (6), the gener- 
alised correlation integral can be defined as: 



C g (r) 



1 

NM 



(10) 



where M is the number of centres, N is the total number of parti- 
cles included in the sample, rii(r) is the same expression defined 
in the Equation (7) and q is called the structure parameter, which 



<(^-<^» 2 > _ (F 2 q ) 



(14) 



If for any qi 7^ q 2 , D qi = D q2 is verified, it is said that the 
distribution is a homogeneous fractal (mono-fractal). For q > 1, 
D q explores the scaling behaviour in high density environments 
(clusters and superclusters) and for q < 1 values, D q explores the 
scaling behaviour in low-density environments, i.e., voids (Sarkar 
et al. 2009). In the event that the distribution of dark matter haloes 
undergo the transition to homogeneity, all values of the fractal di- 
mension must close to the physical space dimension (D q — > 3) and 
the spectrum of lacunarity must tend to zero ($ q — > 0), at the same 
radial distance r. 

Next we will apply the concepts of average mass-radius di- 
mension, generalised dimension and generalised lacunarity to the 
dark matter haloes spatial distribution generated by the Millennium 
dark matter simulation. 



3 MILLENNIUM SIMULATION 

The Millennium Simulation (Springel et al. 2005) is one of the 
most important computational efforts in contemporary cosmology. 
It uses the most accepted parameters inside the standard cosmol- 
ogy, with particles evolving since redshift 2 = 127 until the present 
z — 0. Although the standard cosmological model contains two 
ingredients which have not yet been verified by laboratory experi- 
ments (the components of dark matter and dark energy), the ACDM 
model is almost universally accepted by cosmologists as the best 
description of the present observational data. The spatial geome- 
try is very close to flat and the initial perturbations which should 
have given origin to galaxies are Gaussian, adiabatic, and nearly 
scale-invariant. The cosmological parameters assumed inside this 
computational algorithm are summarised in Table 1. 

Here the density parameters Qj — pi / p CT u are defined from 
the critical density p crit = SH^/SnG. The total present mat- 
ter density composed by dark matter and baryonic matter Q. m — 
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fi m fl\ h erg Ti 

0.25 0.045 0.75 0.73 0.9 1 
Table 1. Cosmological parameters for the Millennium Simulation 
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Figure 1. Millennium Simulation Fractal Analysis: Number of dark matter 
haloes vs Radial distance over 100 centres randomly chosen, in logarithmic 
scale. The initial data dispersion and the slope growth is shown. 

&b + O dm , the dark energy density Q.a with Q. m + = 1, the 
Hubble parameter h = Ho/ (W0kms~ 1 Mpc" 1 ), as the linear- 
theory power spectrum variance in spheres of radius 8 Mpc/h and 
the spectral index n used to describe the density perturbations in 
Fourier space, where n = 1 corresponds to Harrisson-Zeldovich 
scale invariant spectrum. 

The simulation is implemented with the length of the side 
of the simulation cube Lb ox = 500Mpc/h, the number of par- 
ticles N part — 2160 3 « 10 10 , with the resulting particle mass 
m p — 8.6 x 10 8 M Q //i and the gravitational softening length 
rj = 5kpc/h. The Virgo consortium, provides free access to the 
database, which authorised users can entry using the Structured 
Query Language (SQL). It is expected that different researchers 
employ the information to explore how the structures of dark mat- 
ter and galaxies evolve in the standard cosmology. 

We use the information provided by this simulation in order to 
perform a fractal analysis of the dark matter haloes clustering, aim- 
ing to find the way in which the dark matter haloes are grouped as 
a function of the radial distance and to detect the scale of transition 
to homogeneity. 



4 FRACTAL MASS-RADIUS DIMENSION AND 

LACUNARITY OF MILLENNIUM SIMULATION DATA 

The first step in fractal analysis of dark matter haloes clustering, for 
z = from the millennium simulation, is the choice of a random 
sample of centres in the entire simulation (~ 1.5 x 10 7 dark matter 
haloes). We chose randomly one hundred centres with the purpose 
to determine the functional relationship between the logarithm of 
the number of dark matter haloes and the logarithm of the radial 
distance. We want to know if there is a linear function that leads to 
the determination of the fractal correlation dimension D m and the 
pre-factor F as shown in Figure 1 . 

In the case of a mono-fractal set, the expected relationship be- 
tween the logarithm of the number of particles and the logarithm 



of radial distance should be close to a linear behaviour. By observ- 
ing Figure 1 a large scatter in the relationship for distances r < 20 
Mpc/h is evident. Also, it would be inappropriate to assign a single 
slope for the entire curve when a gradual increase of slope with the 
radial distance is observed, suggesting a fractal dimension growing 
with scale. It is therefore necessary to calculate the fractal dimen- 
sion and the pre-factor F in small sections of the curve for each 
centre, in order to determine if there is a radial distance dependency 
of the fractal quantities. 

Dividing the curve in segments using the sliding window tech- 
nique (Martinez & Saar 2002; Rodrigues, Barbosa & Costa 2005) 
and making a linear least-square fit through each set of successive 
points in the log-log plot, it is possible to calculate the slope and 
intercept of the line that approximates the curve section over every 
point for each of the curves generated from each centre, and then 
determine the average over a thousand centres (more representative 
sample) with their respective standard deviation. We also analyse 
the behaviour of this average mass-radius fractal dimension as the 
sample size approaches to the limits of the simulation, as shown in 
Figure 2. 

In Figure 2, the graphs for different depths have a similar 
behaviour; a rapid increase in the fractal dimension to radial dis- 
tances about 20 Mpc/h followed by slower growth until it reaches 
the physical dimension of space to depths beyond the 100 Mpc/h. 
Those graphics with radial distances from the centres exceeding 
180 Mpc/h, where the dark matter haloes are close to the edges of 
the simulation, show that the average mass-radius fractal dimen- 
sion at short distances is larger than the physical dimension of the 
space where the whole fractal is embedded, a physically unaccept- 
able situation. 

The behaviour of the pre-factor F is shown in Figure 3. The 
values of the pre-factor change with scale, starting with high val- 
ues near the centres and diminishing to reach a low constant value 
at distances larger than 100 Mpc/h. Therefore, the average dis- 
tance between dark matter neighbours decreases with radial dis- 
tance reaching a constant value for radial distances larger than 100 
Mpc/h. A change at short distances of the behaviour of the curves 
for depths from the centres larger than 180 Mpc/h is detected. 

This results confirm that there are spurious homogenisation ef- 
fects in the calculations of the fractal dimension. Thus, it is neces- 
sary to calculate the fractal dimension over radial distances smaller 
than 180 Mpc/h, i.e., radial distances where the fractal dimension 
does not exceed 3 ± la, as recommended by Wen et al. (1989). 
In our case the maximum radial distance that satisfies the above 
requirement is 160 Mpc/h. 

To complement the vision of the dark matter haloes fractal 
distribution we calculate the lacunarity up to 160 Mpc/h radial dis- 
tances, taking into account the pre-factor of the sliding window 
least-square fit analysis to the same depth within the simulation. 
The lacunarity of the dark matter haloes distribution, Figure 4, 
shows a maximum value near each chosen centre, followed by os- 
cillatory lacunarity behaviour at radial distances smaller than 60 
Mpc/h. In addition, above 60 Mpc/h we see a smoothly decreas- 
ing lacunarity function with tendency to homogeneity at radial dis- 
tances « 120 Mpc/h where lacunarity values approach to zero. 

The fractal quantities show a strong reliance on the scale. The 
decrease in the lacunarity curve and the scale dependence of the 
average mass-radius fractal dimension with the radial distance sug- 
gest us that the clustering of dark matter haloes do not behave like 
a mono-fractal. Therefore it is necessary to study the dark matter 
haloes distribution from the multi-fractal viewpoint. 



© 2011 RAS, MNRAS 000, 1-13 



Multi-fractal and Lacunarity in Millennium DM haloes 5 



r 



20 40 



Fractal Dimension to 1 40 Mpc/h h 
Physical Space Dimension 



r(Mpc/h) 



Fractal Dimension to 180 Mpc/h y 
Physical Space Dimension 



80 100 
r(Mpc/h) 



120 140 160 180 




Fractal Dimension to 160 Mpc/h i 
Physical Space Dimension 



r(Mpc/h) 



Fractal Dimension to 200 Mpc/h y 
Physical Space Dimension 



20 40 60 



100 120 140 160 180 200 
r(Mpc/h) 



Figure 2. The average of mass-radius fractal dimension vs radial distance for the dark matter haloes distribution in the Millennium simulation. Spurious 
homogenisation effects in the curves around 20 Mpc/h are observed to fractal calculations for depths larger than 180 Mpc/h. ±1(7 error bars are shown. 



5 MULTI-FRACTALITY AND MULTI-LACUNARITY OF 
MILLENNIUM SIMULATION DATA 

In the last section we have shown that the Millennium Simulation 
dark matter haloes average mass-radius fractal dimension have a 
radial distance dependency with a lacunarity descending to zero 
around 120 Mpc/h. Now, the results of the generalised correlation 
integral for twelve different values of q are presented in Figure 5. 
This calculations include overdense regions (g ^ 1) and low den- 
sity regions (q < 1). For q = 1 the numerical limit of the gener- 
alised correlation integral is shown. To avoid the problems of spuri- 
ous homogenisation found in the mass fractal dimension, the multi- 
fractal calculus is limited to depths still far from the borders of the 
simulation. Based on the average mass-radius fractal behaviour, the 
maximum depth for the calculations in the Millennium data is 160 
Mpc/h. 

From the generalised correlation integral, the sliding window 
technique is used to determine the fractal dimension spectrum. We 
find the multi-fractal dimension behaviour as a function of the ra- 
dial distance for the structure parameter values q < 1 as shown 
in the Figure 6. Besides, we calculate the corresponding % 2 test in 
each multi-fractal relation over all distance range where we made 
the calculations; the low % 2 values provide us confidence on the 
used method. The \ 2 average and maximum values for each struc- 
ture parameter are shown in Table 2. 

For structure parameter values q > 1 the radial distance de- 
pendency is presented in Figure 7. In this manner we covered the 
range of the structure parameter — 6 < q < 6 in steps of one. The 
X 2 test was calculated in the same manner as for the multi-fractal 



relations with q < 1. There are low \ 2 values over all radial dis- 
tances, See Table 2. 

In order to detect the transition to homogeneity in terms of the 
fractal dimension, we use the relationship: 

£ 9 -3|/3<e, (15) 

where e is a fiducial value which is calculated from the average of 
the relative dispersion a/D q . In Table 3 we present the shortest 
radial distances for which this condition for the transition to homo- 
geneity is fulfilled. 

This results are confirmed by the relationship between the 
multi-fractal dimension and the spectral parameter g, as shown in 
Figure 8. 

Another possible definition for transition to homogeneity is: 

2[C 2 (r 1 )/V(r 1 ) - C 2 (r 2 )/V(r 2 )] 
[C2(r 1 )/V(r 1 ) + C2(r 2 )/V(r 2 )] ' K ' 

Wi, r 2 > Rh, where is the homogeneity scale, and e is a fidu- 
cial value. If the same value for epsilon as above is used, then the 
transition to homogeneity appears at scales near 20 Mpc/h, which 
do not agree with the scales of transition to homogeneity we found 
from the multi-fractal spectrum or from the minimal values of the 
lacunarity spectrum. Therefore, either this definition for transition 
to homogeneity is less appropriate than the former one or the value 
of epsilon should be much more smaller than 10~ 3 . In the latter 
case, it is not clear what fiducial value should be assigned to e and 
how to estimate it. 

To corroborate this result and because the fractal dimension 
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Figure 3. Fractal Pre-factor F vs radial distance for the dark matter haloes distribution in the Millennium simulation. There is a tendency to achieve a constant 
value with increasing radial distance. ±lcr error bars are shown. 
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Table 2. \ rnean and maximum values obtained for each multi-fractal dimension D q . Top: low density environments —6 ^ q C 0, bottom: high density 
environments 1 ^ q ^ 6. For all the D q analysed the \ 2 values show a peak at scales smaller than 15 Mpc/h and very low, almost constant, values at larger 
scales. 



does not indicate in what form the dark matter haloes set is filling 
the whole space, the analysis is completed with the calculation of 
the spectrum of lacunarity for overdenses regions, as shown below 
in Figure 9. 



6 DISCUSSION 

We confirm the dependence of fractal dimension with the radial dis- 
tance; The most of the previous works reported only a single frac- 
tal dimension value to frame the clustering behaviour of large-scale 
mass (Durrer & Sylos Labini 1998; Sylos Labini 1999; Joyce et al. 
2005). The majority of recent works focus on the radial distance 
dependency of the fractal dimension and the multi-fractal analysis 



of large scale clustering in the Universe, for example: Ramos et al. 
(2002) describe the scaling properties of large-scale structures in 
the Universe and the transition to homogeneity around 500 Mpc/h, 
Seshadri (2005) found that the universe is homogeneous over scales 
larger than about 80 to lOOMpc/h, Bagla, Yadav & Seshadri (2008) 
observed that the fractal dimension makes a rapid transition to val- 
ues close to the physical space dimension at scales between 40 
and 100 Mpc, while Sarkar et al. (2009) determinate that the Mil- 
lennium simulation exhibit a transition to homogeneity at around 
70 Mpc/h. In our case, the analysis of The Millennium Simulation 
confirms fractal behaviour (in the average mass-radius relation) at 
scales smaller than at least 100 Mpc/h with homogeneity transition 
around 110 Mpc/h for the dark matter haloes distribution. 

On the other hand, we observe limitations in the application 
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Figure 4. Mass-radius quantities calculated with the sliding window technique: The average of the mass-radius dimension D m , the pre-factor F, the lacunarity 
<E> and the x 2 test vs the radial distance, for the dark matter haloes distribution. ±lcr error bars are shown. The average < \ 2 >= 7.9 x 10~ 3 , with 
X 2 max = 5.1 x 10- 2 . 
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Figure 5. Generalised Correlation Integral C q (r) for all the values of q studied in this paper for the Millennium simulation dark matter haloes. Over densities 
q 1 left, low densities q < 1 right. All graphics in log-log scale 
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Table 3. Scale of homogeneity transition for every fractal dimension in overdenses environments 1 ^ q ^ 6. The e value is calculated from the relative 
dispersion of the fractal dimension 
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Figure 6. Multi-fractal Dimension Spectrum D q (r) vs radial distance r for low density environments — 6 ^ q ^ from the Millennium dark matter haloes. 
±1ct error bars in the graphics, with 0.99 confidence level in t-student test. 
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Figure 7. Multi-fractal dimension spectrum D q (r) vs radial distance r for high density environments 1 ^ q sg 6 from the Millennium dark matter haloes. 
±1ct error bars in the graphics, with 0.99 confidence level in t-student test. 



of the average mass-radius fractal dimension and the multi-fractal 
dimension spectrum to radial distances larger than 1 80 Mpc/h. We 
detect spurious homogenisation effects from these calculations as 
the radial distance from the centres approaches the edges of the 
simulation. This effect is reported in the literature (Wen et al. 1989; 
Sylos Labini et al. 1996). As stated above, due to this edge effect 
it is necessary to limit the maximum radius to which the fractal di- 
mension can be calculated within the simulation. In our case the 
average mass-radius dimension, the lacunarity, the multi-fractal di- 
mension and the lacunarity spectrum have to be limited to maxi- 
mum radial distances of 160 Mpc/h, where the fractal dimension 
does not exceed the value of 3 ± la. 

The behaviour of mass-radius fractal dimension can be di- 
vided in two regions clearly observed: a fast-growing region from 



each centre up to « 20 Mpc/h in radial distance, followed by a 
slow-growing region to achieve homogeneity at radial distances 
around 100 Mpc/h. Above this distance the search for the scale of 
transition to homogeneity using the lacunarity and the multi-fractal 
treatment should focus. 

The calculation of lacunarity as complement in the determi- 
nation of mass-radius fractal dimension, in addition to confirm 
the process of increasing the homogeneity (Decreasing lacunar- 
ity values), is capable to detect heterogeneity oscillations with lo- 
cal increase of lacunarity. This behaviour indicates alternation be- 
tween voids (lacunarity high values) and concentrations of mat- 
ter (lacunarity low values) and was reported earlier by Provenzale, 
Spiegel & Thieberger (1997), Solis & Tao (1997), Murante et al. 
(1998). Furthermore, the high values of lacunarity at radial dis- 
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Figure 8. Spectrum of the Multi-fractal dimension D q (r) as a function of structure parameter q in the range [—6, 6] for Millennium simulation dark matter 
haloes. Homogeneity transition around 120 Mpc/h is determined. The solid line is an interpolation by the Bezier method. 
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Figure 9. Multi-lacunarity spectram «3> 9 (r) in logscale vs radial distance r for high density environments q > 1 from the Millennium dark matter haloes 
distribution. Oscillatory behaviour modulated by a decreasing function. ±l<x error bars are shown, with 0.99 confidence level in t-student test 



tances smaller than 60 Mpc/h confirm the existence of voids near 
each centre, followed by oscillatory lacunarity behaviour (gaps and 
high concentration of dark matter alternatively). 

The multi-fractal spectrum calculation exhibits less dispersion 
than that found on the average mass-radius fractal dimension and 
its corresponding lacunarity, as seen in the behaviour of the \ 2 test. 
This gives more confidence on the subsequent analyses. The results 
for structure parameter values q < 1, show an excessive increase 
of the fractal dimension for radial distances smaller than 40 Mpc/h, 
reaching values that overcome the physical dimension of space (un- 
physical values), followed by a tendency to homogeneity for voids 
at larger radial distances. In the work of Sarkar et al. (2009), on 
transition to homogeneity in the Millennium Simulation data, it is 
shown that for q = —2 fractal dimension values larger than 3 ap- 



pear not only in the region below 40 Mpc/h, but also a radial dis- 
tances between 70 Mpc/h and 90 Mpc/h. 

For q > 1 we observe a slower growth of the multi-fractal 
dimension as the parameter structure increases, indicating a het- 
erogeneity growth as hierarchical clustering increases. We find the 
possible homogeneity transition located in radial distances between 
100 Mpc/h and 120 Mpc/h . This result has be corroborated by 
evaluating the multi-fractal dimension as a function of structure 
parameter q, where a homogeneity transition around 120 Mpc/h is 
detected. This homogeneity transition scale is larger than those re- 
ported in similar works such as Gaite (2007), Sarkar et al. (2009), 
but lower than the reported by Yadav, Bagla & Khandai (2010). 
Other works on galaxy surveys did not report any homogeneity 
transition (e.g. Sylos Labini & Pietronero (2008) and Verevkin, 
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Bukhmastova & Baryshev (2011)). For this reason we propose to 
use the spectrum of lacunarity as another homogeneity transition 
indicator. 

The lacunarity spectrum presents a similar behaviour to that 
found in the lacunarity calculations based in the average mass- 
radius dimension. An initial increase in the heterogeneity, i.e., great 
presence of voids in regions close to the centres (below 20 Mpc/h) 
followed by a decreasing function which module oscillations with 
lacunarity values alternation. This shows that the clustering of dark 
matter haloes in the Millennium Simulation have low dense regions 
followed by regions of high concentration of dark matter haloes. 
Furthermore, the lacunarity spectrum points to a homogeneity tran- 
sition (minimum lacunarity values) for radial distances between 
100 Mpc/h and 130 Mpc/h, confirming the region where the cross 
to homogeneity was found in the multi-fractal spectrum. 



7 CONCLUSIONS 

In this paper we have studied the distribution of the Millennium 
Simulation dark matter haloes from a multi-fractal viewpoint. Our 
results provide a complementary method on the detection of the 
homogeneity transition scale based on the lacunarity concept. Our 
main conclusions are as follows. 

The dark matter haloes clustering in the ACDM Millennium 
Simulation have a radial distance dependency divided in two re- 
gions. First a region with a fast dimensional growth which starts 
with great variability in regions close to the centres, around 20 
Mpc/h, a symptom of inhomogeneity at short distances. A sec- 
ond region where the clustering of dark matter shows evidence of 
slower dimensional growth arriving to homogeneity with a progres- 
sive decrease in the dispersion. Similar results are reported by Gaite 
(2005, 2007) in his simulation fractal analysis. Therefore, the use 
of multi-fractal analysis appears indispensable because of the com- 
plexity in the clustering of the dark matter haloes. 

In the determination of the fractal mass-radius dimension and 
multi-fractal spectrum it is necessary to consider the spurious ho- 
mogenisation effects introduced as the radial distance approaches 
the edges of the volume to analyse. In our case when it comes near 
the edges, the number of possible centres are being limited (fewer 
centres from which to estimate the relationship); Besides, the mass- 
radius fractal dimension starts to exceed the spatial dimension that 
the clustering of dark matter haloes is embedding, this can lead 
to confusing behaviour over the cross to homogeneity. (Wen et al. 
1989; Martinez & Saar 2002; Gabrielli et al. 2005). For this rea- 
son the maximum radial distance for the fractal analysis should be 
limited in the calculation, according to the sample spatial size. 

For the case of low-density environments, q < 1, it is remark- 
able the coherence between high dimension values for voids below 
r w 20 Mpc/h, the first region of fast growing in the multi-fractal 
dimension spectrum for q > 1, and the lacunarity spectrum for 
the same spatial region. This is an evidence of complexity due to 
self-gravitation at short distances. 

In the calculation of the multi-fractal spectrum, we show the 
homogeneity cross at depths between 100 Mpc/h and 120 Mpc/h, 
without exceeding the physical dimension of space. The cross to 
homogeneity is more precisely located by the relation between the 
multi-fractal dimension and the structure parameter q, which shows 
the cross to homogeneity at 120 Mpc/h radial distance. The spec- 
trum of lacunarity confirms this result, every lacunarity function 
points the homogeneity scale in the same spatial region. 

The lacunarity spectrum for every structure parameter in over- 



denses environments q ^ 1 reveals regions with relative maxima, 
allowing us to detect spatial regions where voids are forming inside 
the clustering of dark matter haloes. The transition of homogene- 
ity for every structure parameter is founded beyond the 100 Mpc/h. 
Insofar that the structure parameter increases, the minimum value 
of lacunarity is located farther from the centre, in the same manner 
that the scale of homogeneity transition is increasing in the multi- 
fractal spectrum. 
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